Static and rotating domain-wall crosses in Bose-Einstein condensates 
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For a Bose-Einstein condensate (BEC) in a two-dimensional (2D) trap, we introduce cross pat- 
terns, which are generated by intersection of two domain walls (DWs) separating immiscible species, 
with opposite signs of the wave functions in each pair of sectors filled by the same species. The 
cross pattern remains stable up to the zero value of the immiscibility parameter |A|, while simpler 
rectilinear (quasi-lD) DWs exist only for values of |A| essentially exceeding those in BEC mixtures 
(two spin states of the same isotope) currently available to the experiment. Both symmetric and 
asymmetric cross configurations are investigated, with equal or different numbers N ly 2 of atoms in 
the two species. In rotating traps, "propellers" (stable revolving crosses) are found too. A full 
stability region for of the crosses and propellers in the system's parameter space is identified, un- 
stable crosses evolving into arrays of vortex-antivortex pairs. Stable rotating rectilinear DWs are 
found too, at larger vlues of jA|. All the patterns produced by the intersection of three or more 
DWs are unstable, rearranging themselves into ones with two DWs. Optical "propellers" are also 
predicted in a twisted nonlinear photonic-crystal fiber carrying two different wavelengths or circular 
polarizations, which can be used for applications to switching and routing. 
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I. INTRODUCTION 

Domain walls (DWs) separating immiscible species are 
generic dynamical structures in mixed Bose-Einstein con- 
densates (BECs) [1]. Originally, the DWs were studied in 
one-dimensional (ID) BEC models, but they can be nat- 
urally extended into the 2D geometry as quasi-lD objects 
[2,3]; in particular, circular DWs, between a less repul- 
sive component in the middle of the trap and a more 
repulsive one forming an outer shell, have been found. 
Recently, complex 2D large-area structures in rotating bi- 
nary BECs were predicted in simulations, including vor- 
tex lattices and sheets in mixtures [4,5]. 

Our aim is to construct genuinely two-dimensional DW 
patterns, in the form of crosses formed by intersection of 
two DWs, in trapped binary BECs. We will consider 
both symmetric and asymmetric crosses, depending on 
the ratio of the numbers of atoms in the two species 
("stoichiometric ratio"). Patterns in a rotating trap will 
be considered too. A rotating DW cross seems like a 
"double- vane propeller" , which may be also symmetric 
or asymmetric. Note that the DWs may move, in the 
general case; however, their motion is not amenable to 
straightforward observation in the ID settings, because 
of the relatively small size of domains available in BEC 
experiments. Thus, the "propeller" offers a unique pos- 
sibility to create permanently moving DWs. We identify 
stability regions for the quiescent and rotating crosses, 
and investigate the evolution of unstable ones. All multi- 
handed crosses, formed by the intersection of more than 
two DWs, are shown to be unstable. We will also briefly 
consider rotating one-dimensional (rectilinear) DWs, i.e., 
"single-vane propellers" (however, we conclude that the 
rectilinear DWs are less relevant than the crosses for ex- 
periments with currently available BEC mixtures). 

An appealing feature of the cross structures is the fea- 



sibility of their experimental realization in binary BECs. 
Experimental creation of two-component mixtures has 
already been reported for different spin states in 87 Rb 
[6,7] and 23 Na [8]. Accordingly, the use of a mixture of 
two spin states of the same isotope seems to be the most 
straightforward way of creating the DW cross configu- 
ration, as the latter can be imprinted onto the BEC by 
optical beams passed through a properly designed phase 
mask. Moreover, the DW cross appears, in our extensive 
computations, to be the most robust structure among the 
various types of domain walls considered. Indeed, the im- 
miscibility condition for the repulsive BEC mixture (the 
one with positive scattering lengths of atomic collisions) 
is 



ana22 - a 12 < 0, 



(1) 



where an, CY22 and a±2 are strengths of the intra- and 
inter-species interactions, respectively [see Eqs. (5) and 
(7) below]. As is shown below, the crosses (both static 
and rotating ones) remain stable up to |A| =0, while 
the rectilinear DW is stable only for 



A < a££ = -0.061. 



(2) 



On the other hand, experimental measurements [6,7] 
yield an extremely small actual value of |A| in the mix- 
ture of two spin states of 87 Rb, 



A 



Rb 



-0.0009, 



(3) 



which does not satisfy the condition (2). For the a mix- 
ture of different spin states in 23 Na, the immiscibility 
parameter extracted from experimental measurements is 
larger, 



A 



Na 



-0.036 



(4) 



1 



[8], but it does not meet the condition (2) either. 

We also note in passing that work is currently in 
progress towards the creation of two-component BECs 
with different atomic species, such as 41 K: 87 Rb [9] and 
7 Li: 133 Cs [10]. In that case, |A| may take values very 
different from those given in Eqs. (3) and (4). Fur- 
thermore, we should note that an additional possibility 
to push the experimental values of A toward the critical 
one mentioned above, is through the use of the, so called, 
Feshbach resonance [11] that can control the strength of 
the inter-atomic interactions. 

The same "propeller" effect may be implemented in 
completely different physical media, namely, photonic 
crystals (PCs) or photonic-crystal fibers (PCFs), with a 
sclf-defocusing nonlinearity and a superimposed twist. In 
PCFs, the twist bends the inner holes, which run along 
the fiber, into helices. In that case, the cross configu- 
ration may be realized as an intersection of stable opti- 
cal DWs separating two orthogonal circular polarizations 
or two different carrier wavelengths (previously, spatial- 
domain optical DWs between different polarizations were 
predicted in driven dissipative cavities [12], but no exper- 
imental setting was proposed to observe their permanent 
motion). The twist-induced rotation of the DW cross 
may take place along the propagation distance, so that 
the entire pattern will seem like a double helix. Twisted 
PCs per se were studied in a stack model [13], and very 
recently, the first twisted PCF was fabricated [14] (it was 
used as a polarization filter). Note that, combining a 
modulation of the air holes density along the radius of 
the PCF core and a properly structured cladding, the 
distribution of the effective refractive index across the 
PCF can be made similar to the parabolic trapping po- 
tential in the BEC model (see below). 

Besides the first possibility to observe permanent mo- 
tion of optical DWs (in the spatial domain), the pro- 
posed configuration suggests a design of a new all-optical 
switch: the turn of the sectors occupied by a signal field, 
with a specific circular polarization or wavelength, just 
by 90° is sufficient for complete switching. The diam- 
eter of the twist-tolerant PCFs, ~ 100 jmm, and their 
length, a few centimeters (which is tantamount to tens 
of rotation periods) [14], make these experiments and 
applications quite feasible. Moreover, the use of the 
co-propagation of several different wavelengths in the 
PC fiber can make stable multi-handed cross and pro- 
peller possible. The latter may be promising for routing 
applications in wavelength-division-multiplexed (WDM) 
telecommunications. Similar configurations with more 
than two species may also be relevant in terms of immis- 
cible multi-component BEC mixtures. The latter issue 
will be considered elsewhere. 

The paper is organized as follows. In Section II we con- 
struct the rectilinear and cross DW configurations and 
study some of their properties. In Section III, we focus 
on the stability of crosses, including rotating ones. The 
above-mentioned "single-vane propeller" (a rotating rec- 
tilinear DW) are also discussed. Section IV summarizes 



our findings and presents our conclusions. 



II. CONSTRUCTION AND PROPERTIES OF 
DOMAIN- WALL CONFIGURATIONS 

The two-component rotating repulsive BEC is de- 
scribed by a system of coupled Gross-Pitaevskii (GP) 
equations [15], 



lh ot 



9jk\tpk\ 



fe=l,2 



1>j, i = l,2, (5) 



where ipj are the wave functions of the two species, nor- 
malized so that Nj — J\ipj\ 2 dr is the respective num- 
ber of atoms. The single-species Hamiltonian in Eq. 
(5) is H = -(?i 2 /2m)V 2 - lolL z + V, where m is the 
atomic mass (assuming a mixture of two spin states of 
the same isotope), ujl is the rotation frequency, and 
L z = ih(xd y — yd x ) is the angular-momentum operator. 
The trapping potential is 



m 



V=-(cj. 



2 2, 22 

r + uj z z 



). 



(6) 



where r 2 = x 2 + y 2 , and the confining frequencies uj r and 
u) z are assumed to obey the condition u) r /u) z e H < 1. 
The intra- and inter-species interactions are character- 
ized by the coefficients gjj — Airh 2 ajj jm and 512 = 321 = 
Airh 2 ai2/m, respectively, where ajk are the correspond- 
ing scattering lengths; as mentioned above, we consider 
the (most typical) case of positive ajk- Then, the condi- 
tion of the immiscibility between the components takes 
the well-known form of Eq. (1). 

Following Refs. [16], effective 2D GP equations can 
be derived from the 3D ones. To this aim, measur- 
ing the coordinates and time in unit s of the harmonic- 
oscillator length and period, i.e., ^jTi/mLo z and l/w z , 
respectively, we seek for solutions to Eqs. (5) as 
i>j(r,z,t) = (2Tr)- 1 / 4 y /hu z /g 11 u j (r, t)$j{z) exp (-i-yjt), 
where &j(z) = 7r -1 / 4 cxp(— z 2 /2) is the ground state of 
the ID harmonic oscillator. Multiplying the resulting 
equations by <J>* and integrating it in z, we arrive at a 
system of 2D equations, 



dt 



2D 



Y a jk\ u k 
fc=l,2 

H 2D = (-l/2)V 2 ± + V(r)-n L L 



1,2, 



(7) 
(8) 



where V 2 ^ is the 2D Laplacian, and 

V(r) = (l/2)Q. 2 r 2 ,L z = i(xd v - yd x ),n L = u L /u z . (9) 

The nonlinearity coefficients in the 2D system (7) 
are an = 1, a\i = 512/511, and a 2 2 = 522/511- 
Then, the numbers of atoms in the species are Nj — 



2 



(4V27T 3 / 2 ) ^h/mL; z Qj, where Q 3 = f \uj\ 2 cPr± are 
the norms of the 2D wave functions. 

The spatial evolution of bimodal optical beams in the 
above-mentioned twisted PCFs is also described by Eqs. 
(7) for the amplitudes Uj of the two components, with 
t replaced by the propagation distance z, and V(r) in 
Eq. (8) is a potential function describing the cladding 
surrounding the PCF proper. In this case, Nj are total 
powers of the two components of the beam, and the non- 
linear coefficients take values close to a i2 = 2au = 2a 22 ; 
according to Eq. (1), the latter values definitely guar- 
antee the "immiscibility" (mutual repulsion of the two 
optical modes). 

Solutions to Eqs. (7) in the form of either a simple 
rectilinear DW (first, without the rotation, VLl = 0) are 
constructed as follows. We start with the Thomas-Fermi 
(TF) configuration for the two components with chemical 
potentials 

( u j)tf = e-^VW^n^] Mi- ( 10 ) 

if — V(r) > 0, and (uj) TF = otherwise (in the PCFs, 
— fij are the propagation constants of the optical modes). 
Then, to construct a rectilinear DW, we effectively re- 
move each component from a half-plane, upon multiply- 
ing {uj) TF by 

/^^[l-t-F/tanhx]. (11) 

The resulting function u\ (u 2 ) is positive in the left 
(right) half of the (x, y) plane. 

To construct a DW cross in a similar fashion, we 
also start from the TF ansatz (10), and then remove 
two quarter-planes in each component, multiplying the 
expression (10) by (l/2)[f J (x)f J (y) + f 1 (-x)f J (-y)}, 
where the factors /i,2 are defined as per Eq. (11). As 
a result, the wave functions u\ and u 2 have support, re- 
spectively, practically only in the quadrants (sectors) 1 
and 3 (i.e., xy > 0), and 2 and 4 (i.e., xy < 0). Next, we 
multiply the entire configuration by (— l)- 7 tanhec, hence 
the functions u\ and u 2 are, respectively, positive (neg- 
ative) in the sectors 1 (3) and 2 (4). In other words, we 
introduce the phase shift of 7r between the two quadrants 
filled by each species (in other words, we aim to construct 
a DW cross which, simultaneously, is a "latent dark soli- 
ton" in each component) . Without these phase shifts, the 
cross will be obviously unstable against splitting into a 
set of two approximately parallel quasi-rectilinear DWs, 
as a pair of the sectors filled by the same species without 
the phase shift between them will tend to merge into a 
single stripe-like pattern. 

After the prototype rectilinear or cross DW configura- 
tion was constructed as described above, Eqs. (7) were 
numerically integrated in imaginary time, to let the sys- 
tem relax into a stable stationary state closest to the pro- 
totype one (which is a well-known technical ruse). For 
the DW-cross case, the relaxation in imaginary time did 



not essentially alter the prototype pattern, always con- 
verging to a numerically exact stable cross pattern, pro- 
vided that the immiscibility condition (1) was met. This 
outcome of the imaginary-time integration took place for 
arbitrarily small values of the immiscibility coefficient 
|A|, including the case when A was set precisely equal to 
0. With A > 0, the integration always showed that two 
species would completely mix into a uniform state. 

In the case of the rectilinear DW, the imaginary-time 
integration converged to a stable quasi-lD DW pattern 
only if the condition (2) was satisfied, which is more re- 
strictive than the general immiscibility criterion (1). For 
smaller values of |A|, stable rectilinear DWs could never 
be found; instead, the two species would completely mix 
up, even if A was negative. Thus, the DW-cross pattern, 
although being more complex than its simplest rectilin- 
ear counterpart, is a much more robust one, and has a 
much better chance to be created in experiments with the 
currently available weakly immiscible spin-state mixtures 
that actually do not satisfy the condition (2), see Eqs. 
(3) and (4). The extra robustness of the DW-cross struc- 
ture may be explained by the additional tt phase shifts 
lent to its both components, as described above. 

As an example, in Fig.l we show the result of the re- 
laxation for each type of the prototype pattern, with the 
value of A taken as for the spin-state mixture in 87 Rb, see 
Eq. (3). This was realized by setting an — 1, C122 = 0.94, 
and «i2 = 0.97 in Eq. (7). As is seen, the rectilinear DW 
pattern cannot be formed indeed, preferring to mix itself 
up into a uniform state [panels (a) and (b)]. On the con- 
trary to that, the initial (prototype) DW-cross pattern 
readily relaxes into a configuration of exactly the same 
type. 

In the example shown in Fig. 1, the normalized 
magnetic-trap strength is O = 0.05 [see Eq. (9)], and 
both chemical potentials arc set equal to 1. In terms of 
the real- world parameters, this choice corresponds to a 
mixture of two spin states in the Rb condensate in a 
disk-shaped trap with 

lu t . = 2ir x 6 Hz, uj z = 2ir x 120 Hz (12) 

[see Eq. (6)], the initial TF radius and numbers of atoms 
in each species being 

R = 34^m, N = 3.6 x 10 3 . (13) 

In what follows below, we will display examples for 
larger |A|, which may be unrealistically large directly 
(but may be attainable through the use of the Fesh- 
bach resonance [11] mentioned above). This allows us 
to generate patterns that are much sharper and easier 
to understand, while qualitatively they are completely 
tantamount to those found at smaller (more directly re- 
alizable) values of the immiscibility parameter. So, we 
will fix 

an = 1, a 22 = 1.01, ai 2 = 1.52, (14) 



3 



which corresponds to A = —1.3. 
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FIG. 1. Gray-scale plots showing the distribution of the 
density k*i| 2 of one species (the density \u 2 \ 2 is complemen- 
tary to |wi| 2 ) in the initial configuration (a) corresponding to 
the prototype rectilinear DW, and in the corresponding final 
counterpart (b), which is generated by the numerical integra- 
tion of the GP equations (7) in imaginary time, in the case 
of weak immiscibility, with A = —9 x 10~ 4 (corresponding 
to an = 1, Q 2 2 = 0.94, a 12 = 0.97) and fi = 0.05. The 
panels (c) and (d) display the same for the initial and final 
configurations in the case of the DW cross. 
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FIG. 2. Gray-scale plots of the densities of the two species, 
|tti| 2 and |«2 1 2 (left and right panels), in the symmetric (a, 
b) and asymmetric (c, d) DW-cross patterns. The param- 
eters are an = 1, 022 = 1.01, ai 2 = 1.52, and fi = 0.05. 
In the symmetric configuration, jUi = fj,2 = 1, and in the 
asymmetric one = 1.1025, fj,2 = 1, and N1/N2 = 1.2, 
61 = 6-k/U, 62 = 5-k/U. The estimate (in physical units) for 
the Thomas-Fermi diameter for both configurations is 68/xm. 



III. STABILITY AND ROTATION OF THE 
DOMAIN- WALL PATTERNS 



A. Static domain-wall crosses 



Before proceeding further, it is important to note that, 
alongside the symmetric DW-cross configurations like the 
one shown in Fig. 1(d), asymmetric ones are possible 
too, with different numbers of atoms in the two species, 
Qi 7^ Qi- Indeed, the TF approximation (10) yields 



* - (4) - i tffl2 ) fl2 ' (15) 



where R is the TF radius of the state; consequently 
an 7^ a 22 and/or fix ^ [i 2 lead to Q\ ^ Q 2 , and, 
as a result, to asymmetric crosses. We have checked 
that the numerical results always obey a natural rela- 
tion, Ni/N 2 = Q1/Q2 = O1/82, where 6j is the intrinsic 
angle of the j-th sector, with 9\ + 62 = 7r. Examples of 
the symmetric and asymmetric DW crosses are shown, 
respectively, in the top and bottom panels of Fig. 2. 



Dynamical stability of the symmetric and asymmetric 
DW crosses is a crucially important issue. We have found 
that, for a fixed trap's strength 0, the stability strongly 
depends on the stoichiometry ratio Ni/N 2 , more sym- 
metric configurations being more robust. To show this, 
we fix the nonlincarity coefficients ctjj as in Eqs. (14), 
we set fi2 = 1 and, then, for a given value of O in the 
interval (0,0.25), we vary to induce variation of the 
ratio N1/N2 = Q1/Q2, as per Eq. (15). Finally, we sim- 
ulate Eqs. (7) in real time (up to t = 1000), to test the 
stability of the configuration. The resulting stability do- 
main in the (Ni/N 2 , ft) parametric plane is displayed in 
Fig. 3(a). For example, at ft = 0.05, which corresponds 
to the above-mentioned values (12) and (13) of the phys- 
ical parameters, Fig. 3(a) shows that stable DW crosses 
exist in the interval 0.65 < N\/N 2 < 1.65. We note that, 
alternatively, it is possible to represent the domain of 
stability of the static DW crosses in the (^1,^2) plane, 
for a fixed value of ft (then, each point of the domain 
would correspond to a given value of Ni/N 2 ). This can 
be done, e.g., upon rescaling the variables in Eqs. (7)-(8) 
so that = 1 (see also a relevant discussion in Section 
III.B); in that case the domain of stability would be a 
triangular-like zone enclosing the line \i\ = 
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Region III: Instability 
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FIG. 3. Stability regions for the static (a) and rotating (b) 
DW crosses. The values of the nonlinear coefficients are fixed 
as in Eqs. (14). In (b), solid and dashed lines, respectively, 
show borders of the stable-rotation region for symmetric and 
asymmetric crosses (in the latter case, Ni/N 2 = 1.5). 
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FIG. 4. Evolution of the density |ui| 2 of the species with 
the larger number of atoms in an unstable DW cross, for the 
same parameters as in Fig. 2, but with u\ — 1.8025. In this 
case, N-i./N 2 = 2.5 and 9i = 5tt/7, 6 2 = 2tt/7. The snap- 
shots (a)-(f) correspond to t = 0, t = 100, t = 310, t = 316, 
t = 1000 and t = 2200, respectively (the physical time unit is 
1.33 ms). The density \u 2 \ 2 is complementary to |ui| 2 , as in 



Figs. 2(c, d). 



Evolution of unstable DW crosses is also an issue of in- 
terest. A typical example of the instability development 
(in the component with the larger number of atoms) is 
displayed in Fig. 4 for hi = 1.8025 and [i 2 = 1, which 
corresponds to N1/N2 = 2.5. At an initial stage, the 
cross quickly rearranges into a quasi-lD object, which 
is actually a dark soliton in the species with the larger 
number of atoms, coupled to a bright soliton in the other 
species (this object resembles structures considered in 
Ref. [3]). The latter configuration is itself subject to a 
snaking instability, which is a known feature of quasi- ID 
dark solitons in BECs filling 2D areas [18] (we stress that, 
unlike the dark solitons, the rectilinear DWs are not sub- 
ject to this instability). The snaking instability initiates 
break-up of the dark stripe into four vortex-antivortex 
pairs in the first species, with the second species collect- 
ing itself into spots coupled to the pairs. Finally, the 
four vortex-antivortex pairs annihilate into two. The lat- 
ter configuration persists for long times, and it may be 
related to serpentine-shaped vortex sheets reported in re- 
cent 2D simulations of large-size BECs [4] . It would be 
desirable to analyze this instability by means of a finite 
mode stability analysis in the form of [19]. However, in 
the present setting, the modulus of the wavefunction is 
not radially symmetric (nor does the DW configuration 
bear a complex phase structure), hence it is not directly 
amenable to such an approach. Understanding the origin 
of such an instability would consitute an interesting topic 
for future studies. 



B. Domain- wall crosses in the rotating trap 

For the rotating trap, with fir, 7^ in Eqs. (7), we 
have found that both the symmetric and asymmetric DW 
crosses revolve stably in a certain frequency interval, 

(^) min <«L<(^) max . (16) 

If fi/, < (^L) m ; n , the cross docs not rotate at all, while 
at fi/, > (^i) max it decays (see below). These results 
are summarized in Fig. 3(b). In region I, the DW cross 
remains quiescent, in region III it gets destroyed, while 
the stable rotation occurs in region II. Note that there are 
minimum values of the trap strength which are necessary 
for the rotation: fi m i n = 0.012 and 0.018 for the symmet- 
ric and asymmetric "double- vane propellers" . Those val- 
ues correspond to the edge points on the solid and dashed 
lines, respectively, in Fig. 3(b), the rotation frequencies 
at these points being fi/, = 0.0034 and 0.010. 

Note that the horizontal axis in Fig. 3(b) is extended 
up to a relatively large value, fi = 0.25 (recall that it was 
assumed above that the magnetic-trap's strength fi is a 
small parameter, in order to derive the 2D GP equations 
from the underlying 3D system). According to experi- 
ence accumulated in applications of asymptotic methods 
to various models (including the GP equation), the value 
0.25 is small enough to produce reliable results. 

Typical values of the physical parameters admitting 
the stable rotation of the DW cross can again be esti- 
mated for a mixture of two spin states in the 87 Rb con- 
densate. For instance, taking the trapping frequencies 
Lu r = 2n x 7 Hz and lu z = 2n x 70 Hz, the TF radius 
R = 22 ^m, and 10 3 atoms in each species, the rota- 
tion frequency ul = 2tt x 0.84 Hz definitely falls within 
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the stable-rotation interval (16). Examples of a stably 
rotating propeller with N\/N 2 — 1.2, and of its self- 
destruction in the case of Ql > (^i) max (for Ni = N 2 ) 
are shown in Figs. 5 and 6, respectively. 
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FIG. 5. A stable rotating assymetric domain-wall cross 
( "double- vane propeller"). The four snapshots (a), (b), (c), 
and (d) show the distribution of the density |mi| 2 at consecu- 
tive time moments - respectively, t = 0, T/6, T/3, and T/2, 
which cover a half of the rotation period, T/2 = ~ 52.3 

(~ 70 ms, in physical units). The parameters are as in Fig.2(c, 
d) but with Q = 0.1 and n L = 0.06. 




FIG. 6. An example of the evolution of an unstable sym- 
metric "double- vane propeller" (rotating domain- wall cross). 
The panels in the left and right columns display the density 
distributions \ui\ 2 and \u2\ 2 , respectively, at t = [frames 
(a) and (&)], t = 10 [frames (c) and (d)], and t = 20 [frames 
(e) and (/)]. The parameters are as in Fig. 2(a,b), but with 
fi = 0.1 and Q. L = 0.08. 

As it was mentioned above, a dynamical property of 
the DW crosses which is crucially important for their rel- 
evance to the currently available experimental settings, 
that have a very small value of the immiscibility param- 
eter |A| [see Eq. (1)], is the fact that the crosses remain 
stable for all the negative values of A, up to A = 0. This 
robustness carries over to the rotating crosses, as illus- 



trated by Fig. 
the case of A = 



7, that shows a stable rotating cross in 
: 0. 
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FIG. 7. A stable rotating symmetric domain-wall cross 
found exactly at the immiscibility border, A = (correspond- 
ing to an = «22 = CK12 = 1) and = 0.1. The arrangement 
of the figure is the same as in Fig. 5. 



The results can alternatively be summarized upon con- 
sidering a rescaling of variables in Eqs. (7)-(8). In partic- 
ular, upon measuring the time, spatial variables and nor- 
malized wavefunctions in units of SI -1 / 2 and fi 1 / 2 
respectively, we obtain a system of two equations similar 
to Eq. (7), but with the Hamiltonian being given by, 



(17) 



Then, we may present the domain of stability of the rotat- 
ing DW crosses in the parameter plane (Ni/N 2 ,£Il/£1) 
as follows: We set /l*2 = 1 and, then, for a given value 
of fii (which sets the value of Ni/N 2 ), we vary the ratio 
in the interval (0, 1). Then, we numerically inte- 
grate the system of equations for Uj, for long times, to 
test the stability of the rotating DW cross. The results, 
are shown in Fig. 8, where the stable-rotation region is 
delineated in the full parameter plane: In region I, the 
DW cross remains quiescent, in region III it is destroyed, 
while the stable rotation occurs in region II. Note that for 
values of Ni/N 2 outside the shown interval, (0.65, 1.55), 
the DW crosses are unstable. Thus, for instance, at a 
typical value of the ratio between the rotation frequency 
and the trap strength, fi^/fi = 0.1, the stoichiometry ra- 
tio must belong to the interval 0.65 < Ni/N 2 < 1.55. As 
mentioned above, the smallest value of Q^/Q necessary 
for the stable rotation (which is 0.0019) corresponds to 
the symmetric case, Ni/N 2 = 1. 
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FIG. 8. The region of the stable rotation of the domain-wall 
crosses, shown in the parameter plane (N1/N2, fli/fl). In re- 
gion I, the DW cross is quiescent, in region III it is destroyed, 
while in region II it is rotating in a stable manner. Note that 
for values of N1/N2 outside the interval (0.65, 1.55) the DW 
crosses are unstable. 



C. More general cross and propeller patterns 

All the patterns produced by the intersection of more 
than two DWs were found to be unstable, rearrang- 
ing themselves into the fundamental crosses considered 
above. A typical example of that, with four intersecting 
DWs and Ni/N 2 = 1, is displayed in Fig. 9. The same 
instability of the higher-order crosses takes place in the 
rotating trap. 



y 



y 



y 




-20 x 20 -20 x 20 

FIG. 9. Instability of a higher-order domain- wall cross (left 
and right panels show distributions of the densities |iti| 2 and 
|ii2 1 2 ) - The frames (a) and (b) display the initial configuration 
at t = 0. Further frames show the evolving configurations: 
(c), (d) at t = 50, and at (c), (f) at t = 600. The parameters 
are as in Fig. 2(a, b). 

Two-dimensional patterns similar to the symmetric 



DW crosses can be formed in single- component BECs by 
two linear dark solitons intersecting under the right an- 
gle (provided that the dark soliton itself may be stable). 
We have checked that such single-component patterns are 
always unstable. 

The crosses in BECs (but not in the above-mentioned 
photonic-crystal media) can be extended into the 3D 
case, as patterns formed by the intersection of planar 
DWs. Moreover, torque can be applied to such a struc- 
ture, making it look like a double helix. Further, these 
patterns may revolve together with the trap. Results for 
the cross-shaped 3D patterns will be reported elsewhere. 



D. Rotating rectilinear domain walls 

As was explained above, stable rectilinear DWs are less 
relevant the currently (experimentally) tractable small 
values of the immiscibility parameter |A|. Nevertheless, 
for completeness of the analysis, and also for the sake of 
possible future experiments in BEC mixtures of different 
atomic species, where |A| may be much larger than in 
the presently available mixtures of different spin states 
of the same isotope, it makes sense to briefly consider a 
possibility of stable rotation of the DWs of this type too. 
Numerical simulations show that, generally, the rectilin- 
ear DW withstands rotation, within a certain interval of 
the angular velocities [cf. Eq. (16)], if it is stable in 
the static trap. Here, we do not aim to produce compre- 
hensive results for the stability domain of the rotating 
rectilinear DWs. Instead, in Fig. 10 we display an exam- 
ple of a stably rotating DW of this type. This example 
has a special purport, as it shows the rotating rectilinear 
DW found at the minimum value of the immiscibility pa- 
rameters at which it may be stable, |A| = 0.061. Stable 
rotating rectilinear DWs found in the generic case are 
quite similar to the one shown in Fig. 10. 




FIG. 10. Rotation of the rectilinear domain wall at the 
smallest value of the immiscibility parameter, |A| = 0.061 
(corresponding to an = 022 = 1, Qi2 = 1.03), and £2 = 0.1 at 
which it is stable. The arrangement of the figure is the same 
as in Fig. 7. 
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IV. CONCLUSIONS 

Our analysis predicts the existence of symmetric and 
asymmetric domain-wall (DW) crosses in immiscible 
mixtures of repulsive BECs in disk-shaped traps. The 
stable crosses are distinguished by the phase shift of tt 
between sectors filled by each component. The same pat- 
terns may also be carried by bimodal light beams in pho- 
tonic crystal fibers (PCFs). The crosses are stable in a 
wide range of parameters, including a certain interval of 
the stoichiometry ratio Ni/N 2l the case N\ = N 2 being 
the most robust one. Adding the trap's rotation (or twist 
of the PCF, in the optical model), we have shown that 
the DW crosses can rotate, provided that the driving an- 
gular velocity (or the twist pitch, in the PCF) is limited 
from below and from above (if it is too small, the cross 
does not rotate, and if it is too large, the cross decays). 
In the PCF, the rotation takes place not in time, but in 
the spatial domain, i.e., along the propagation distance, 
hence the optical cross actually looks like a double he- 
lix. The rotating optical crosses have potential for the 
use in switching and routing applications. Full stabil- 
ity domains for the quiescent and rotating crosses were 
identified, corresponding to values of physical parame- 
ters accessible in current experiments (both for BECs 
and PCFs). 

The crosses (including the rotating ones) persist ex- 
actly up to the zero value of the immiscibility parameter 
| A|. This aspect of the robustness of the DW crosses is 
crucially important because, in the BEC mixtures of dif- 
ferent spin states in Rb and 23 Na, currently available to 
the experiment, the actual value of | A| is very small. We 
have demonstrated that, in the cases of practical interest, 
the simplest rectilinear (quasi-lD) DWs do not exist (in 
a stable form), while the crosses remain entirely stable. 
For the completeness' sake, we have also demonstrated 
a possibility of stable rotation of the rectilinear DWs at 
larger values of |A|. 

The patterns produced by the intersection of more 
than two DWs were observed to be unstable. They rear- 
range themselves into the fundamental crosses. 

Finally, it is worthy to mention that, at angular veloc- 
ities of the rotating trap essentially exceeding the range 
of values dealt with in this paper, one may expect the ex- 
istence of rotating DWs separating not merely areas with 
the quasi-uniform distribution of the densities of immis- 
cible species, but rather areas filled with triangular vor- 
tex lattices (each species supporting its own lattice). In 
single-component BECs, such regular lattices, composed 
of a large number of vortices, have been recently stud- 
ied in detail in direct experiments [20]. In fact, creation 
of the DW separating the vortex lattices in immiscible 
components is the most straightforward way for experi- 
mental observation a rotating DW pattern in BECs. This 
possibility will be considered in detail elsewhere. 
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